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APPLICATION OF ERTS-1 DATA TO ANALYSIS OF AGRICULTURAL CROPS 
AND FORESTS IN MICHIGAN 

Gene R. Saflr, Wayne L. Myeri, Michigan Agricultural Experiment Station, Michigan 
State University, East Lansing, Michigan; William A. Malila, Jamei P. Morgemtern, 
Environmental Research Institute of Michigan, Ann Arbor, Michigan 


ABSTRACT 

The results reported are based on analysis of ERTS Fri me 1033- 
15580 collected over southwestern Lower Michigan on August 25, 1872, 
Major agricultural crops such as corn and soybeans were approaching 
maturity at this data and forest canopies were dense. 

Extensive ground truth information was gathered by detailed field 
study of test strips. This detailed information was supplemented over 
larger areas by Interpretation of RB-57 and C-47 photography and MSS 
imagery. The S. D. A. - A. S. C. S. also cooperated by providing infor- 
mation on crops from their records. 

Recognition processing of ERTS-1 MSS data was carried out on a 
digital computer.- Fields and forest stands were selected as training 
sets and test areas. Aerial imagery was essential for locating the 
positione of these selected areas on ERTS digital tapes. 

The recognition process was successful for each type of vegetation 
which had a dense green canopy such as forests, corn, and soybeans. 
Bare soil was also recognizable as a category. However recognition of 
species was difficult in senesclng or senescent vegetation. Since the 
accuracy of recognition depends cm stage of growth, optimum times for 
coliftcting data will vary from one crop to the next. 

Accurate computer recognition of crops from satellite data will 
be useful in operational surveys as the first stage in a multistage 
sampling process. 
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Introduction 


Michigan State University (MSU) in cooperation with the Envi- 
ronmental Research Institute of Michigan (ERIM) began a program in 
the summer of 1972 to tea', the usefulness of ERTS-1 satellite data for 
monitoring and managing crops and forests in Michigan. Specifically, 
the objectives included: (1) verif' cation that major agricultural crops 
and forest types can be identify . ■ rom ERTS-1 data; (2) development, 
application, and testing for accuracy of multispectral techniques for 
crop and forest acreage estimation in Michigan; (3) correlation of 
' ariations in signatures from space with ground truth data. 

In addition io the scientists directly involved in analysis of the 
data, a team of cooperat s has been assembled to evaluate the 
operational utility of the results which emanate from the project. 

This team includes members of state, federal, and local agricul- 
tural and natural resource agencies. 

Tata analysis to date has been confined to ERTS frame E-1033- 
lf iu (August 25, 1972) as a result of the inclement weather condi- 
tions which prevailed throughout most of M'chigan's 1972 growing 
season following the launch of ERTS-1. 

Ground Truth Information 


Direct field observation and 35-mm photography were the main 
sources of ground truth information for the analysis of agricultural 
crops. Specifically, biological parameters such as plant height, row 
direction and width, percent ground cover, stage of growth, corn 
tassel color, and disease incidence were estimated and recorded 
for numerous selected fields in the test area. Since forest cover 
changes less rapidly than agricultural crops, the primary source of 
information for forests was photointerpretation of RB-57 and C-47 
underflight imagery. The photo 1 rTorpretive work was supplemented 
by collection of data on the ground as necessary. The RB-57 and C-47 
imagery was also extremely -iseAil for analysis of agricultural crops. 

In addition to the field and underfiight data, cooperators in both the 
United States Department of Agriculture, Agricultural Stabilization 
and Conservation Service (U. S. D, A. -A. S. 0 C- i and the Forest Service 
have contributed to the pool of ground truti mfor mutton. The A. S. C. S. 
efforts produced a set of annotated copies oi r .urged airphotos showing 
the location and nature of vegetation types on the holdings of landowners 
who subscribe to A. S. programs. 
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M. S. S. Digital Analysis - Methods 


Digital tape data for frame 1033-15580, were screened for quality 
by preliminary processing on the ERIM digital computers. They were 
found to exhibit the same problem present in a set of tapes for the same 
frame received by ERIM under another contract. The ’ roblem is that 
reproduced signals from one of six detector elements which generate 
the MSS data in ERTS band 6 (0. 7-0. 8 pro) are faulty. Thus, anomalous 
data are present for band 6 in every sixth line of data; otherwise, the 
data appear to be satisfactory. This problem complicates signature 
extraction and data analysis. In particular, recognition processing 
for the work described here was restricted to three channels. 

The priina.., st sites (in Eaton Co. , Michigan) were located 
within the digital data, and lLie-printer gray maps were produced for 
all ERTS bands. The gray maps for ERTS bal'd 5 were used to locate 
selected training and test plots of known ground cover. The RB-57 and 
C-47 underflight imagery was essential for correctly locating these 
plots, which wei e then designated by line and point number to the com- 
puter for extraction of signal statistics. In the selection of training 
sets, care was taken to avoid boundary points. Fifty-eight lots wer. 
designated and ERTS signal statistics were extracted for eight types of 
ground cover. These statistics were subjected to cluster analysis, 
and the results were used to select several plots for combination 
to form recognition signatures. The plots which were not used di- 
rectly for specifying signatures became "test" sets for evaluating the 
accuracy of recognition. Eighteen additional test plots were then sel- 
ected and included in the analysis. 

Recognition maps were produced for an intensive test area in 
Eaton Co. , Michigan. Recognition runs were based on the three good 
ERTS channels using several different sets of parameters. Firs*, 
twelve recognition signatures were used and map3 were produced with 
different rejection threshold levels. That is, each observation was 
classified as belonging to one of the recognition signatures and then 
tested to see if it was unlikely enough to be rejected and categorized 
as belonging to none of the classes considered. Next, seven recogni- 
tion signatures were used; the seven recognition signatures included 
combinations of the pairs of signatures used for several classes in 
the twelve ignature runs. 

MSS Digital Analysis - Results 

Recognition results were analyzed for the 76 identified plots. The 
overall results of the first-look analysis of recognition are summarized 
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in Table 1 for five cover classes (corn, soybeans, trees, senescent 
vegetation, and soils). As noted earlier, only three ERTS channels 
were used (4, 5 and 7). The values in Table 1 represent averages of 
percentages computed separately for each plot analyzed. The overall 
average percentage of correct classification (for test sets) is over 
83%. ^he average percentage error is 10%, with 16% being Type 1 
(i. e. , missed classification, including not classified) errors and 4% 
being Type II (i. e. , incorrect classification) errors. If "not classified" 
points are excluded from the computation, the overall average is 
nearly 85% correct. 

Recognition percentages are high for those vegetation classes that 
had mature and uniform canopies at the time the data were collected 
(Aug. 25th). Corn, soybeans, and trees(forest) met this criterion, 
and were classified accurately. The class of senescent or senescing 
vegetation included observations from field beans, wheat stubble, and 
grass. These canopies were characterized by non-uniform distribu- 
tions of dead a’-d dying vegetation along with patches of more healthy 
vegetation. For example, field beans had matured and had begun 
senescing, while soybeans and corn were more vigorous. Also, wheat 
stubble fields were dry and brown except for some that had been 
seeded to alfalfa or red clover; the latter fields had patches of green 
growth among the stubble. The wide variability within these vegeta- 
tion types at this time of year makes it difficult to classify them 
accurately. Alfalfa is a crop that is harvested repeatedly at irregular 
intervals throughout the growing season, and plots of it can appear 
very different, depending on their conditions at the time of observation. 
One vigorous alfalfa was included initially and accurately recognized. 

A lack of test plots, for which the exact condition at the time of the 
ERTS-1 pass is known, caused us to omit alfalfa as a ciass from the 
reported analysis. Bare soil was distinctive and accurately recognized. 

Thus, the first -look analysis for computer recognition within 
boundaries of selected plots shows a good capability for differentiating 
each type of vegetation that has a dense green canopy, with bare soil 
also being recognizable as a category. The next step in the analysis of 
computer recognition is a more critical evaluation of accuracy by 
cover type for all resolution elements in selected portions of the frame. 

E lement-by-Element Analysis for Forest Cover 

Figure 1 is a portion of the gray map for ERTS band 5 in Chester 
and Roxand Townships of Eaton County, Michigan with major roads 
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delineated. An RB-57 color infrared photo was used to transfer the 
locp*.ions of the forests to the gray map, and the elements that fall with- 
in -ne forest area are shown by heavy dots in Figure 1. Figure 1 is 
conr ervative in that most doubtful border elements were not designated 
as "forest". Figure 2 is a computer recognition map for the same 
area. Heavy dots have beei superimposed on the elements which were 
correctly classified as "trees" (forest) . The forest c’ements which 
were not recognized as such by the computer (Type 1 error) are in- 
dicated with a heavy square having a white center. Type II errors 
(incorrectly classified as "trees") are indicated by triangles. The 
Type I error for forests on this portion of the frame is approximately 
40%. An examination of Figure 2 shows that most of these errors take 
place in border elements. For the most part, these border elements 
were classified as corn. The remaining Type I errors are mostly 
accounted for by areas in which the forest canopy is sparse. The 
Type H errors are only about 3%. 

Since the original "trees" training sets were located in *'.e center 
of dense woodlots, the misclassification of sparsely stocked areas is 
not too surprising. An examination of the likelihood for the mis class- 
ified elements showed a very low probability of classification under the 
"trees" signature. Use of separate training sets and subresolution 
element analysis are being investigated as possible means of improving 
recognition in sparse forests. The current classification would give a 
reasonable estimate of the acreage that is suitable for woodlot manage- 
ment, but would give an underestimate for total acreage of forest. 

Summary 

Computer analysis of ERTS-1 data provided good recognition of 
vegetation classes that had mature and uniform canopies at the time 
when the data were collected. Bare sou was also recognized accurately. 
Classification was extremely difficult for senescent vegetation which 
was characterized by non-uniform distribution of dead and dying vege- 
tation along with patches of more healthy vegetation. Since the accuracy 
of classification depends on the stage of growth, optimum times for 
collecting data will vary from one crop to the next. However, the 
optimum for recognizing each crop is yet to be determined. This 
'ears further study, especially for field beans since Michigan is the 
leading producer of this crop in the United States. 





TABLE I. SLMMARY OF RECOGNITION RESULTS ON A PLOT-BY- 
PLOT BASIS FOR 76 PLOTS, ERTS FRAME 1033-15580, 

3 CHANNELS (ERTS 6 EXCLUDED), 0.001 PROBABILITY 
OF REJECTION 


A. 


Average Percentage of Class' Plots Assigned to 

Listed Recognition Signature 

No. No. Senesc 


Class 

Plots 

Points 

Corn 

Soy 

Alf 


Tree 

Bl an 

Grass 

Soil 

Corn 

21 

481 

84. 27 

0. 55 

0. 

13 

9. 85 

3.85 

1. 35 

0 

Soy 

10 

115 

1.00 

89. 40 

2. 

30 

2. 59 

2. 61 

0 

0 

Trees 

12 

358 

11. 00 

3. 80 

0 


84. 50 

0. 20 

a 50 

0 

Seresc 

16 

306 

16. 30 

6. 55 

7. 

15 

0 

54. 30 

8.23 

6. 53 

Soils 

4 

56 

0 

0 

0 


0 

0 

0 

97. 62 

TOTALS 76 

1416 










Summary of Percentages (Averaged Over Plots) 


Class 

No. 

Plots 

No. 

Points 

Not 

Clas'd 

Correctly Incorrect. 
Assigned Assigned 
To Class To Class 

Average 

Error 

Correct 
Excluding 
Not Clas'd 

Corn 

21 

481 

0 

84. 27 

7. 29 

11. 51 

84. 27 

Soy 

10 

115 

2. 10 

89. 40 

. 50 

6. 55 

91. 31 

Trees 

12 

358 

0 

84. 50 

3.66 

9. 58 

84.66 

Senesc 

19 

340 

0.94 

62. 53 

2.59 

20. 03 

63. 12 

Soils 

14 

122 

2. 38 

97. 62 

2.00 

2. 19 

100. 00 

Averaged Over 
Five Classes 

1.08 

83. 66 

0.61 

9. 97 

84. 67 
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Figure 1. Channel 5 gray map for portions of Chester and Roxand 
Townships in Eaton County, Michigan showing actual 
locations of woodlots (•) and roads 


179 


* X. H -4h « II/* ••)••• 












****** »iii^ <* o BH0 * * * * 0 * *00000 Qrf I *•*»--» 
■ 0889* «•* *0000080* • *•>%•*%•*• 0OOB0%#O0BH * *0* *• 0000* •*• BS* -*■ * 
©§0*. •• 1 1 ••• I •• • i **H**OO00****S**-“« 

•O8*****JL%OO0#*******«6* j | *0IMM*i*»lYy*8***#«^»*» | *•♦»--- 

l*i j j •••■■OOOe* • « * * !*»****✓* * • ** j 

•aft*** I 1 1 • 0*«8* 9809*98* ***•••••••••*•«*** ••'Wrff** *• *** **88 - 

S8>***»«*m I t»r* *»»* laaooa 

**8« 1 1 ********** *89*08* *(>•••••••••••8* *** • I ••* ••*•#** lofMI 

• ••••I | ()•••••* BB800B0000OOOO0000OOO*****8 •>*«• *B0«B 

****** I I i.«*B*e»6800BBMMBBBBB»«BMBB***e* = -*** I I** I'h^OMB 
t«SB*************HIIIHS**INHaHI*8***ei«- • *II(***TVU>B 

•••Bee***ee*ieB 8 bo«aB«B****BMMM 38 e****afB — *•*♦ h n***>« 

***080****9BOOO8B0OOO POO * j 

)- -*l 


OOOOOOOOOOOOSMOfOOB •••••## 0000000* -* i 

•000000000000** 1 000****** ••* I 


0«****8«8***** I * *8* ***•**•♦•** *800* *** I ** 
B*****0B***BB000-bB« **80**** I ***BB8*** *•* , 
B **** *88* •*****■-■ 0----* ***** *9B • «**jjm#iT' 
****«**0******B«--- »-«• »** • | 

• •••• |*0**t*«**« — »*»«* 


f 

I 


# - 
U- 


80*1111*** 

eoooooooooooiioooiioe******ooooooooo -*! mooe** 

■ * *-•• — 

. . (a a 

- * 

* **0* | *- — 0 
1 1 *• | * **08**0039 
•B»***** *8890 

******* | **000 

• •• l*!*********** • * 6 1 0«8B** * * * 1 1 ! | ******** |***** | *080 

III ****888«juk*"%'0'V***B**8****B9B****** ij • ****8*** ****** I *10i 

j i 1* *ft#*^8**«6**8BB — BBBBbB8bB****** i**00000f0#08B 

X Wo^8****8*»O80B***B«a = —00008*8* ****•#.*• 0000 0000 000*****9e 

. * * 0600098 *98000 0B 0« 68*068888* * * * * *V oOOOOfiOOOOOO* * • 1 80 

**#**88*8** • *0000808 — «-»0*** , •••8******^ai000i0i000009***B8 
*•**6000 1 i • • 8BOOOBO8880* ■ -88 * « B8B&** * * * * *^*000000000* • ***** 90 

• •••• *00* ***008000* *000*9008* **0000**»***’*B0800* 8 **********0 

?• ♦ • • 888800880008 • • 

iOOO9***>*****30088B00B*** 

•***eO00BO*8«**00000000000i00** *** **8808* *i* I ***80080080* *90 
00000 8f ****0808088090** 

• • • 080000000s*© 

’ll****. |***«i 1 *•0000000*080 *60 1 • •••**B*****#|****0**B*00B808 



0BV*809OOOOB***08********B*** • • •••0B|h^**MTir**O*****« ■••0000 
••l*860iOO***** B ••«•*•**«•*•• | | *t>* *0* • %* ••••••• *8001 

•00i»»o0iie88e«**e* • •* *BajkB«f^r** ** * • *Bb • i * *800Sti* « bh*» * * * *••• 

•08|*690000B8«8**f^» B *OOT*#*#g8 B 00**08*** i B8HB88| **0B **•#*## 
,«#0808*1O*n}^“ #*##*»#8**» 8B«**0W0*****0***1 1 .*R0B*- *|**« 
•••••***••• • B 0 **B******B*H**«***H 869 %’ HB»»B **»00 
[•••••• *6* •••••B***60* 1 8PB* 1 1 ***8** 1*080880%* *8*08* *880 

• **00**t*« ••••••• ***e i i ••*•••••• h* i ••••••• m •••••«• ** o 

J*898BBB8B08**OO08**BB8OO0******* *H6000** •#*%••# III' ‘ 

• • V*00000800***OOB***BOOOiO0****0**oeOOOOOO**Vl**** I 

••i«**V*oooooooefi**o«o***eoooo»o»*»**«**ooooooo***^*** • 

• •«•••% *08000000 **0000* • • 80000000* ***OoOOOOOOOO*0*\***** l . 

eeBBBHBk* • *OO0MO »* 00000 * *80000000* 00000000000 *# *00 ••*••* 00 

0* 0080*1*0000000 00 OOOOOBBOOOOOOOO •• 08> >'■ 080 •••••• *v* ••* *800 

060eeeBi**eBootoeo«oo«oooooMOOO8e*888 800**09* 0^*0009000 

• •8B000WJ *6000090* *OOOOOO60OOOOO0B**0B* “•*“0000008008% 0*ri 3000 
0B080H*e0*0*08****«OOOOO**O«OOO8BO**e0B “"0*080090 88 H0000000 


**e 

:« 

|*H 


Figure 2. Computer recognition map showing correctly and 

incorrectly classified elements associated with forest 
cover; • indicates correctly classified forest areas, 
a indicates forest areas mlBclassified, and % indicates 
non-forest areas classified as forests. 
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